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Abstract 

Making use of an ansatz for the eigenfunctions, we obtain an exact closed 
form solution to the non-relativistic Schrodinger equation with the anharmonic 
potential, V(r) = ar 2 + 6r~ 4 + cr~ 6 in two dimensions, where the parameters of 
the potential a, b, c satisfy some constraints. 
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1. Introduction 



The exact solutions to the fundamental dynamical equations play crucial roles in 
physics. It is well-known that the exact solutions to the Schrodinger equation are 
possible only for the several potentials and some approximation methods are frequently 
applied to arrive at the solutions. On the other hand, in recent years, the higher 
order anharmonic potentials have attracted much more attention to physicists and 
mathematicians [1-3]. Interest in these anharmonic oscillator-like interactions stems 
from the fact that, in many case, the study of the relevant Schrodinger equation, for 
example in atomic and molecular physics, provides us with insight into the physical 
problem in question. 

Recall that in the three-dimensional spaces, rough speaking, there are two main 
methods to be used to deal with the anharmonic potentials V(r) = ar 2 + br~ 4 + cr~ 6 . 
One [4, 5] is based on an ansatz for the eigenf unctions to obtain an exact solution with 
this potential. This method undoubtedly provides an exact solution for the ground 
state but sometimes with some constraints on the parameters of the potential. The 
other [6, 7] is relied on a Laurent series ansatz for the eigenfunctions, which converts 
the Schrodinger equation into a difference equation and then the continued fraction 
solutions are defined. This method, however, does not give any constrains for the 
parameters of the potential. 

The reasons why we write this paper are as follows. On the one hand, with the 
advent of growth technique for the realization of the semiconductor quantum wells, the 
quantum mechanics of low-dimensional systems has become a major research field. Al- 
most all of the computational technique developed for the three-dimensional problems 
has already been extended to lower dimensions. On the other hand, the study of the 
potential V(r) = ar 2 + br~ A + cr~ 6 in two dimensions has never been appeared in the 
literature. We now attempt to research it in two dimensions. 

This paper is organized as follows. In Sec. 2, we study the ground state of the 
Schrodinger equation with this potential using an ansatz for the eigenfunctions. The 
first excited state will be discussed by the same way in Sec. 3. The some constraints 
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on the parameters of the potential a, b, c are given in Sees. 2 and 3. The figures for the 
unnormalized radial functions are plotted in the last section. 



2. The ground states 



Throughout this paper the natural unit h — 1 and /j, — 1/2 are employed. Consider 
the Schrodinger equation with a potential V(r) that depends only on the distance r 
from the origin 

fid d 1 d 2 \ 

\r or or r l dip 2 J 

where the potential 



V(r) = ar 2 + br 4 + cr 6 , a > 0, c> 0. 



(2) 



Owing to the symmetry of the potential, let 

^(r, v ) = r - x l 2 R m (r)e ±m ^ m = 0, 1, 2, . . . , 
where the radial wave function R m {r) satisfies the radial equation 



d 2 R m {r 
dr 2 



+ 



E - V(r) - 



m — 1/4 



R m {r) = 0, 



(3) 



(4a) 



where m and E denote the angular momentum and energy, respectively. For the 
solution of Eq. (4a), we make an ansatz[4:, 5] for the ground state 



R m o{r) = exp[p m0 (r)], 



(5) 



where 

Pmo(r 

After calculating, we arrive at the following equation 



-ar 2 + -j3r 2 + nlnr. 



d 2 R m0 (r) 
dr 2 



d 2 Pmo(r) ( dp rn0 (r) 



dr 2 



dr 



RUr) = 0. 



Compare Eq. (4b) with Eq. (4a) and obtain the following set of equations 



(6a) 



(46) 



a 2 = a, j3 2 = c, 



(7a) 
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k 2 - K-2af3 = m 2 - 1/4, (76) 

3(3 - 2(3k = b, (7c) 

E=-(2n + l)a. (7d) 

It is easy to obtain the values of parameters for p m o(r) from the Eqs. (7a) and (7b) 
written as 



a = ±-\/a; (3 = ±v^cj K — 77 =t V m2 + 2^/ac. (8) 

In order to retain the well-behaved solution at the origin and at infinity, we choose 
positive sign in k and negative signs in a and (3. According to these choices, the Eq. 
(7c) leads to the following constraint on the parameters of the potential, 

(6 + 2v^) 2 -4c(m 2 + 2^) = 0. (9) 

The eigenvalue E, however, will be given from Eq. (7d) as 

E = M 4 + 7?)- (10) 

Now, the corresponding eigenfunctions Eq. (5) can be read as 



Rmo(r) = N r K exp 



-\(V^r 2 + yfir 2 ) 



where N is the normalized constant. Here and hereafter k = 1/2 + \Jm? + 2^/ac. 



3. The first excited states 

With the same spirit, we make an ansatz for the eigenfunctions corresponding to 
the first excited state in the potential (2) as follows 

R m i(r) = f m (r) exp[p ml (r)], (12) 



with f m (r) given by 



and p m i(r) given by 



f m (r) = ai + a 2 r 2 + a 3 r 2 , (13) 



Pmi(r) = ]j a i r ' 2 + 2 + K i lnr - (66) 



For short, it is readily to see from Eqs. (12) and (13) that the radial function R m i(r) 
has the following relation 



R m i{r)" 



Pmi{r) +(p m i(r)) +1 j-^ 



R m i(r)=0, (4c) 



where the prime denotes the derivative of the radial function with respect to the variable 
r. Calculating Eq. (4c) carefully and comparing it with Eq. (4a), we obtain 



a 2 [E - Va(2Ki + 5)] = 0, a 3 [6 - y^(2 Kl - 7)] = 0, 



(14a) 



ai[£ - Va(2«i + 1)] = a 2 [m 2 - 1/4 + 2^fac 3«i - 2], (146) 

ai[m 2 - 1/4 + 2v/ac - k\ + m] = a 2 [6 - >/c(2«i + 1)] + a 3 [-E + Va(2«i - 3)], (14c) 

ai[b - y/c(2m - 3)] = -a 3 [m 2 - 1/4 + 2^ - k\ + 5«i - 6], (14c?) 

aii = ±>/a, A = ±-\/c. (14e) 

Hence, if the angular momentum m of the first excited state is the same as that of the 
ground state, we obtain from Eq. (14) 

6 + 7^ 



E l = v^(5 + 2ki), 



(15a) 

(156) 
(15c) 
(15d) 



2^ ' 

ai = 0, a 2 = -v/a") ^3 = —y/c, 
f3\ = —y/c, ol\ = — \fa, 
b = —6y/c, 

where the constants both a± and f3\ are chosen negative signs in order to hold well- 
behaved nature of the solution at r — * and r — > oo. Equation (15d) is another 
constrain on the parameters of the potential. 

At last, the eigenvalue Ei and eigenfunctions R m \(r) for the first excited state with 
the potential (2) may be read from Eqs. (15a) and (12) as follows 

6 ' 



£i = Vo 12 + 



(16) 



R m i(r) = N^r 2 + a 3 r- 2 )r K1 exp 



--(v^r 2 + v^r- 2 ) 



(17) 



where Ni is the normalized constant for the first excited state and K\ is given by Eq. 
(15a). 



As a matter of fact, the normalized constants N and Aq can be calculated in 
principle from the normalized relation 



Considering the values of the parameters of the potential, we fix them as follows. 
The value of parameter a is first fixed, for example a = 1.0, the values of the parameter 
c and b are determined by the constraints Eq. (9) and Eq. (15d) for m — 0. By 
this way, the corresponding parameters turn out to a = 1.0, c = 4, b = — 12, k = 
— 1.5, K\ = 0.5,02 = 1,03 = — 2. The ground state and the first excited state energies 
corresponding to these values are obtained as E = —2 and E 1 = 6, respectively. 
Actually, when we study the properties of the ground state and the first excited state, as 
we know, the unnormalized radial wave functions will not affect the main features of the 
wave functions. We have plotted the unnormalized radial wave functions R l , (i = 0, 1) 
in fig. 1 and fig. 2 for the ground state and the first excited states, respectively. 
Comparing them with the figures for the ground state and the first excited state in 
three dimensions, respectively, it is easy to find that they are different from each 
other. The reason is that the parameters of the potential b, c are not same as those 
in three dimensions, which origins from the different constraints on the parameters of 
the potential, even if the parameter a is same both in two dimensions and in three 
dimensions. 

To summarize, we discuss the ground state and the first excited state for the 
Schrodinger equation with the potential V(r) = ar 2 + br~ A + cr~ 6 using a simpler 
ansatz for the eigenfunctions and simultaneously two constrains for the parameters of 
the potential are arrived at from the compared equations, which then results in the va- 
riety for the energy eigenvalue and eigenfunctions with the varieties of the parameters 
of the potential. This simple and intuitive method is easy to be generalized. The other 
studies to the sextic potential and the octic potential as well as the inverse potential 
in two dimensions are in progress. 




(18) 
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